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Abstract 

We demonstrate how the uncertainty of parameter point estimates can be 
assessed in a maximum hkehhood framework in order to prevent overfit- 
ting and erroneous detection of time-inhomogeneity. The class of models we 
consider are regular vine (R-vine) copula models, for which we describe a 
new algorithm for the exact computation of the score function and observed 
information. R-vine copulas constitute a flexible class of dependence mod- 
els which are constructed hierarchically from bivariate copulas as building 
blocks only, and our algorithm exploits the hierarchical nature for subse- 
quent computation of log-likelihood derivatives. Results obtained using the 
proposed methods are discussed in the context of the asymptotic efficiency 
of different estimation methods for R-vine based models. In a substantial 
application to a dataset of exchange rates, we obtain clear indications for 
time-inhomogeneous dependence between some currency pairs. 

Keywords: copula, exchange rates, rolling-window analysis, R-vine, 
standard errors, time- variation 



1. Introduction 

The last years have seen the rise of dependence models for complex mul- 
tivariate data and the renewed experience during the recent (2008) finan- 
cial crisis that well- accepted economic paradigms may loose their validity 
from one day to another. This has created significant inte rest in th e stud y 
of time-homogeneity of multivariate dependence (s ee e.g^ 



Manner and Reznikoval (l201l[ ) , IChoUete et al.l ( 120091 ) and lGarcia and Tsafack 



Pelletiei (2006) 



(120111 )). However, while it is a standard exercise in multivariate statis- 
tics to compute the uncertainty incorporated in parameter point estimates 
for classes like the multivariate normal distribution, this is often not pos- 
sible for the more complex models which are required to accurately cap- 
ture the dependence in real-world multivariate data. A particularly suc- 
cessful class of dependence models are regular vine (R-vine) co pula model s 
which have been int r oduced in a s ubsequent series of papers by Jod (19961) 



Bedford and Cookd (l200lL 120021) iKurowicka and Cookd (I2006f ). lAas et al. 
(120091 ) and iDifimann et al.l ( 120111 ). They have been applied to model de- 
pendence in variou s areas including agricultural science a nd electricity loads 
Smith et al.l (120101 ) . exchange rates (ICzado et al.l (120121). IStober and Czado 
(1201 ll )). order books and headache data lPanagiotelis et al.l ( 120121 ). In general, 
an R-vine model is constituted by a set of bivariate copulas corresponding to 
conditional distribu t ions d etermined by a sequence of linked trees. Following 
Bedford and Cookd (l200l[ ). the trees (Ti, . . . ,Td-i) forming an R-vine tree 
sequence V are required to fulfill the following properties: 

1. Ti is a tree with nodes A^^i = {1, . . . ,d} and edges Ei. 

2. For i > 2, Tj is a tree with nodes Ni = and edges Ei. 

3. If two nodes in Tj+i are joint by an edge, the corresponding edges in Tj 
must share a common node (proximity condition). 

An 8-dimens ional example is given in Figure [1] Following the notation of 
Czadol (l2010h with a set of bivariate copula densities 

B = {cj(e),fc(e)|D(e)|e G -Ej, 1 < 2 < c/ - 1 } Corresponding to cdges j (c) , /c(e) | D (c) 
in Ei, foT 1 < i < d — 1 the density of a d-dimensional R-vine distribution 
with marginal densities /fc, 1 < /c < ci is given by 

fl,...,d{Xi, . . . , Xd) 
d d-1 

= n '^j{e)Me)\D(e){Fj(^e)\D(e){Xj(^e)\XD(e)),F^e)\D(e){Xk(e)\XD{e)))^ 

i=l e£Ei 

(1) 



1=1 



Here, Xi:){e) is the subvector of x determined by the set of indices in D{e), 
which is called conditioning set while the indices j(e) and k{e) form the 
conditioned set. If all marginal densities are uniform, the corresponding 
distribution is called an R-vine copula. Often it will be convenient to split the 



2 



specification of bivariate parametric copula families in B from the respective 
parameters which are stored in an additional vector 0. In this case, an R-vine 
copula is specified in terms of RV = (V, B, 0). 

Despite the wide range of applications of R-vine copula based models in prac- 
tice, there is a surprising scarcity in the literature considering the uncertainty 
in point estimates of the copula parameters. It is well known that maximum 
likelihood estimates On will be strongly consistent and asymptotically normal 
under regularity conditions on the bivariate building blocks, i.e. 



On — 



N{0, Idp) as n — cxD. 



(2) 



Here, 6 is the true p-dimensional parameter vector, n is the number of (i.i.d.) 
observations, Idp the p x p identity matrix, and 



1(6) = -Ee 



dOidOj 



d_ 



mx) 



09] 



mx) 



(3) 

denotes the Fisher Information Matrix with l{6\x) being the log-likelihood 
of parameter 6 for one observation x. Given the fact that full maximum 
likelihood inference is numerically difficult when non-uniform marginal dis- 
tributions are involved , also two-step procedures have been developed. In 
particular, iJoe and Xul (119961 ) employ the probability integral transform us- 
ing parametric marginal distributions which are fitted in a first step to obtain 
uniform (copula) data o n which the copula i s estimated using ML in a second 



step. As an alternative, iGenest et al.l (|1995[ ) proposed to use non-parametric 
rank transformations in the first step. While these methods are compu- 
tationally rnore tr actable, they are asymptotically less efficient. Following 
Hobaek Hafla (120 111 ), we can decompose the asymptotic covariance matrix for 
the estimates of dependence parameters in a marginal part and a depen- 
dence part: 

•^0,two—step -y/- dependence _j_ -^margins 

where Y'^'^^'^^'^'^^'^'^''^^ = I{0)^^ for ML estimation and the second part is 
zero only when there is no uncertainty about the margins. Further, also the 
estimation of copula parameters can be performed in a tree by tree fashion. 
Denote the log-likelihood arising from parameters in tree i by 



kiOi\x) :-- 



E 

eeEi 



log {cj(^e),k{e)\D{e){Fj{e)\D{e){Xj(e)\xD{e)),Fk{e)\D{e){Xk{e)\xD{e))\0e)) , 
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where the set of components of 6i is {6e\e G -Ej}. In particular Oi is estimated 
by maximizing li{6i\x) and the obtained estimates Oi are used to calculate 
the arguments of the copula functions in l2{02\x). Now, we maximize l2{02\x) 
in 62 to obta in 62 and proceed until all parameters are estimated (for an 
algorithm see IStober and Czadd ( 120111 ) ) . This implies that h also implicitly 
depends on the parameters 6j for j < i through its arguments, i.e. li{6i\x) = 
li{6i, Oi^i, . . . , 0i\x). The parameter estimates obtained from this sequential 
procedure have asymptotical covariance 



^dependence, seq. ^~^1^Q (jT"^^) 



(4) 



where J'q involves second derivatives of the R-vin e log-likelihood furi ction 
and ICe involves elements of the score function, see iHobaek Hafli ( 120111 ). To 
be more precise, Kq and J^e are defined as 



//C 



\ 



6»,1,1 

0^ 
0^ 



\ 



^e,d-2,d' 
0^ 



(5) 



e,d-l,d-l I 



( 0,1,1 



Je 



9,d' 

e,d- 



2,1 
1,1 



04- 



2,d-2 
14-2 



with Ke 



Je 



-E 



E 



a/,(0,,...,0i|x) 



d9, 
.,0i|X) 



9/j(0j,...,0i|X) 



90, 



T 



\ 



J 04-14-1 J 



and 



(6) 



dOidO, 



.d-1. 



While this asymptotic theory is well known, it is almost never applied in 
practice since the estimation of the asymptotic covariance matrix will in- 
volve the Hessian matrix, i.e. the second derivatives of the R-vine likelihood 
function. For these derivatives, no analytical expressions have been available 
creating a gap between theoretical knowledge about estimation errors and 
practical applicability which we will fill in this paper. 

The remainder is structured as follows: Section [2] shows how the (log-) like- 
lihood of a general R-vine model can be calculated efficiently. Based on the 
illustrated algorithm. Sections [3] and H] consider the computation of the first 
and second derivatives of the R-vine copula log-likelihood with respect to the 
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parameters or the score function and observed information, respectively. In 
Section [S] we set to answering the headhne question of this paper the devel- 
oped algorithms to a dataset of exchange rates and Section [6] concludes by 
wrapping up our results. 



2. Computation of the R-vine likelihood 



In order to calculate the (log-) likelihood function of an R-vine model, we 
must develop an algorithmic way to evaluate the copula terms in the decom- 
position ([1]) with respect to the appropriate arguments. Here, the R-vine 
structure with the proximity condition implies that for each edge e the term 



F(Xj(e)|a3D(e)) 



_ 9Cj(^e)J'{e)\D{e)\j'{e)iF{Xj(^e)\XD(e)\j'{e)),F{Xf(^e)\XD(e)\j'{^^^^ 
dF{Xf(e)\XD{e)\j'{e)) 
hj(e),j'{e)\Die)\j'ie)iF{Xj(e)\XD{e)\j'ie)),F{Xj>(e)\XDie)\j'{e))), 



(7) 



and similarly F{xk(e)\xD{e)) can be computed. This means, that there is 
an index j'(e) G -D(e), such that the copula Cj(e)j'(e)|D(e)\j'(e) is in B. For 
this expression, the assumption that the copula Cj(^e),j'{e)\D{e)\j'{e) does not 
depend on the values XD{e)\j'{e) is made. This is called simplifying assump- 
tion, as it will simplify further computations as well as model selection. For 
which classes of multiv ariate distributions this assumption is or is not ap- 



plicable is discussed in IStober et al.l (l2012l ). To further ease notation, we 



will assume that all copulas under investigation are symmetric in their ar- 
guments, such that we do not have to differentiate between ^j{e),j'(e)|D{e)\j'(e) 
and hj'(^e),j{e)\D{e)\j'{e)- While this is valid for most common parametric copula 
families, we can easily drop this assumption later. 

In order to perform computations for a general R-vine copula model, it is con- 

venien t to use a rnatrix n otat ion which has be e n intr oduced by lMorales-Napoles et al 
( I2OIOI ) , iDifemannI ( I2OIOI ) and iDifimann et al.l ( 1201 ll ) . It stores the edges of an 
R-vine tree sequence in the following way: Consider for example the R-vine 
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in Figure [T] which can be described in matrix notation as follows: 



M : 



■m2,i 
m3,i 

"^6,1 



"12,2 
m.3,2 
"14,2 
"15,2 
"16,2 
"17,2 
"18,2 



"13,3 
"14,3 
"15,3 
"16,3 
"17,3 
"18,3 



"14,4 
"15,4 
"16,4 
"17,4 
"18,4 



"15,5 
"16,5 
"17,5 
"18,5 



"16,6 
"17,6 
"18,6 



"17,7 
"18,7 



"18,8 



7 7 

2 2 

3 3 
6 4 

1 
5 



4 
1 



1 1 



(8) 

For example the edge index 74 1 156 is stored by m2,2, '^5,2 given m6,2, ^^7,2 and 
m8,2- Accordingly, we can store the copula families B and the corresponding 
parameters 6. 



4,m6,5|"l7,5,"l8,l 

^4,m7,5|m8,6 



^3,1717,6 |"18, 6 

^3,m8,6 



[72, m8 



,7 / 



^2,4|1,3 

^3,4|1 ^2,3|1 



71 2 



As an illustration for how the R-vine matrix is derived from the pictured 
tree sequence in Figure [1] and vice versa, let us consider the third column of 
the matrix. Here we have 6 on the diagonal, and 2 as a second entry. The 
set of remaining entries below 2 is {3,4, 1,5}. This corresponds to the edge 
2, 6|1, 3, 4, 5 in T5 of Figure [H Similarly, the edge 3, 6|1, 4, 5 corresponds to 
the third entry 3 in the third column, 1, 6|4, 5 to the fourth entry, etc. Note 
that the diagonal of M is sorted in descending order which can always be 
achieved by reordering the node labels. From now on, we will assume that 
all matrices are "normalized" in this way as this allows to simplify notation. 
Therefore we have irii^i = d — i + 1. For further shortening and clarity of index 
labels, we will illustrate the notations in the example of an 8-dimensional 
vine. 

Applying the matrix notation, the (log-)likelihood of an R-vine model is 
computed by first storing all required conditional distribution functions eval- 
uated at a (i-dimensional vector of observations {ui, . . . , Ud) in two matrices. 
In particular, we calculate 

\ 

U2 Ui ) 

(9) 



V 



direct 



TO7,5 ; "ms,5 



F{u^\u 



U3 



6 




^ ^ 7,811,2,3,4.5,6 ^ ^ 

T^ Q2,8|l,3,4.5.(r) Q2, 7|1, 3, 4. 5. (T) 

Figure 1: An R-vine tree sequence in 8 dimensions with edge indices corre- 
sponding to the pair-copulas in an R-vine copula model. 
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V 



indirect 



F{Ums7\^2) 



-'"Is, 6 



-'"Is,? 



(10) 



Note that, for each pair-copula term in ([T]), the corresponding terms of V'^^^^'^^ 
and can be easily determined. When being able to evaluate 



)) 



we do also obtain 



F(u4|u 



"16,5 5 ' ^mg 



.5) 



= (92^^)4, me, 5|m7, 5, "18,5 (-^(''^4|^m7,5 ) 'Um8,5 



•^"16, 5 I """17,5 ' ^ms^hl )l 



)), 



where (diC) and (820) denote partial derivatives with respect to the first 
and second argument, respectively, c.f. Equation ([7]). With all such condi- 
tional distribution functions being available, the copula terms in ([T]) corre- 
sponding to the next tree T4 can be evaluated. This sequential calcul ation 
is performed in Algor ithm 12. 11 which was developed in DifimannI ( 2010 ) and 
Difimann et al.l ( 120111 ). Following the notation in ([7]), we write h{-, ■\B''^\e''^') 
for the conditional distribution function corresponding to a parametric fam- 
ily S*^'* with parameter 6''^'*, where B^'^ and 6^'^ denote the {k, i)th element 
of the matrices B and 6, respectively. Exempli gratia. 



F{u4\u2,U3,Ui) = h{F{u4\u3,Ui),F{u2\u3,Ui)\B4^2\3,l,0i. 



2|3,lJ 



direct „,indirect\ t26,5 q6,5 



6,5 



6,6 



^b,5^^6,5^)_ 



The only question which is left to solve for the computation of the (log-) like- 
lihood is whether the arguments in each step (i.e. F{u4^\u3,Ui), F{u2\u3,Ui) 
in the example) have to be picked from the matrix V'^^^'^'^^ or y^"-d-^'^^<^^-_ For 
this, we exploit the descending order of the diagonal of M. From the struc- 
ture of V'^^^'^'^'^, we see that the first argument of the copula term with family 
S^'* and parameter 6''^'* is stored as the {k,i)t\i element v"^"^^^ of To 
locate the second entry, let us denote M = {rhk^i\i = 1, . . . ,d; k = i, . . . ,d), 
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where rhk^i '■= max{mfc^j, . . . , j} for all i = 1, . . . , (i and k = i, . . . ,d. The 
second argument, which is "n^fc+i.i, • • • ■>'^d,i) niust be in column 

[d — ihk^i + l) of or ^^^^^^^rf j-^y ^j^^ ordering of variables. If ihk^i = mk,i, 

the conditioning variable Um^ ^ has the biggest index and thus the entry we 
are looking for must be in y"'*''^'^*. Similarly, if rhkA > rrik^i, the variable with 
the biggest index is in the conditioning set and we must choose from 

Example 2.1 (Selection of arguments for C2,4|i,3^ 

As an example for how this procedure selects the correct arguments for copula 
terms in the regular vine let us consider the copula C2,4|i,3 in our example 
distribution. The corresponding parameter ^2,4|i,3 is stored as 6^'^ , thus we 
are in the case where i = 5 and k = 6. Since rfig^s = max{m6,5, ^7,5, "^8,5} = 
max{2, 3, 1} = 3 and ffiQ^^ = 3 > 2 = mg^s we select as second argument 
the entry v^^^Y-^^, +1) = y^^^irect _ i?(^i2|ui, U3). Together with vf^^^^ = 
F{u4\ui,Us) which we have already selected, this is the required argument. 

The corresponding algorithm to compute the log-likelihood of an R-vine spec- 
ification for a single observation u = {ui, . . . , uf) is given in Algorithm 12.11 



3. Computation of the score function 

In this section we develop an algorithm to calculate the derivatives of the 
R-vine log-likelihood with respect to copula parameters and thus the score 
function of the model. Throughout the remainder, we will assume that all 
occurring copula densities are continuously differentiable with respect to their 
arguments and parameters. Further, we assume that the copula parameters 
are all in R, the extension to two or higher dimensional parameter spaces is 
straightforward but makes the notation unnecessarily complex. 
To determine the log-likelihood derivatives, we will again exploit the hier- 
archical structure of the R-vine copula model and proceed similarly as for 
the likelihood calculation. The first challenge which we must overcome to 
develop an algorithm for the score function is to determine which of the cop- 
ula terms in Expression ([1]) depend on which parameter directly or indirectly 
through one of their arguments. Following the steps of the log-likelihood 
computation and exploiting the structure of the R-vine structure matrix M, 
this is decided in Algorithm 13.11 

Knowing how a specific copula term depends on a given parameter, we can 
proceed with calculating the corresponding derivatives. Before we explain 
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Algorithm 2.1 Log-likelihood of an R-vine specification. 

Require: d-dimensional R-vine specification in matrix form, i.e., M, B, 0, 
set of observations {ui, . . . , Ud). 
Set L = Q. 

Let V^^^"* = (t'fcf ^> = l,---,d;k = i,...,d). 

Let yi^direct ^ (^^^^^f rect |^ ^ 1, . . . , d; k ^ i, . . . , d) . 

Set (<rs<rs'--,<r) = 

Let M — {mk.i\i = 1, . . . ,d;k = i, . . . ,d) where rhk^i — 
max{mfc..j, .... m,j,i} for alH = 1, . . . , and k = i, . . . ,d. 
for i = d — 1, . . . ,1 do {Iteration over the columns of M} 
for k = d, . . . ,i + 1 do {Iteration over the rows of M} 
Set zi = vfj''^ 
if rhk,i = mk,i then 

Set ^2 = 
else 

Set Z2 = vV^fr^^ 
end if 

Set L = L + c{zi,Z2\B'''\9''''). 

Set f^!!!";* = Z2\B'''\ e^'') and ^^i^* = h{z2, zi\B^^\ 0^^'). 
end for 
end for 
return L 



the derivatives in detail let us start with an example where two of the three 
possible cases of dependence on a given parameter are illustrated. 

Example 3.1 (3- dim) 

Let xi ~ Fi,X2 ~ -^2,3^3 ~ F3 andul — F{xi),U2 — F{x2), us — F{x3), then 
the joint density can be decomposed as 

/i23(a;i,a;2,a;3) = fi{xi)f2{x2)h{x'i) ■ 01,2(^1, M2|^i,2) ■ 02,3(^2, M3 16*2,3) 

• Cl,3|2(/il,2(M2,Ml|6'l,2), /i2,3(M3,M2 16*2,3) 1^1,312) 
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Algorithm 3.1 Determine copula terms which depend on a specific param- 
eter^ 

The input of the algorithm is a d-dimensional R-vine matrix M with elements 
{''^id)id=K--4 number k and column number i corresponding 

to the position of the parameter of interest in the corresponding parame- 
ter matrix d. The output will be a matrix C (with elements (Qj)ij=i,...,d) 
of zeros and ones, a one indicating that the copula term corresponding to 
this position in the matrix will depend on the parameter under considera- 
tion. 

1: Set g (m^^i, m^.i, . . . , rud^i) 

2: Set := l,j — l,...,d 

3: for a = i, ... ,1 do 

4: for b = k, . . . ,a + 1 do 

5: Set h := {ma,a, rUb^a, "^6+l,a, • • • , md,a) 

6: if #(^n/i)==#^then 

7: Set Cb,a := 1 

8: end if 

9: end for 
10: end for 
11: return C 



The first derivatives of In /123 with respect to the copula parameters are 

a(ln/i23(3;i,3:2,X3)) ^ d0,^Ci^2{Ul,U2\9i^2) 
<96'i,2 Ci,2(mi,M2|6'i,2) 
, 5lCl,3|2(/il,2(M2,Ml|^l,2),/i2,3(M3,M2|^2,3) 1^1,312) o , / m \ 



de 



2,3 



^62,3^2,3(^*2, ^i3 


6*2,3) 


C2,z{u2,Uz 


6*2,3) 



, 92Cl,3|2(/^l,2(^i2,^il|^l,2),/i2,3(^i3,^i2|6'2,3)|^l,3|2) , . 1^ ^ 

H 77 — 7 1^ — r-r — 7 — rr^ ^ " Ce2,3'^2,3lM2, Ms t^2,3j 

Cl,3|2 (aIi,2 (^i2 , ^il 1 6'l,2) , Al2,3 (% , "2 I fc'2,3) I t'l,3|2 ) 

9(ln /i23(a:;i, 3:^2, x^)) de^^^^^Ci^3\2{hi,2{u2, mi|^i,2), /i2,3(M3, '"21^2,3)1^1,312) 



96'i,3|2 Ci,3|2(/il,2(li2,1il|6'l,2),/i2,3(li3,1i2|6'2,3) 16*1,312) 



The first case which occurs in our example is that the copula densities Ci^2 
and C2,3 depend on their respective parameters directly. For a general term 
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involving a copula cu,v\z with parameter 9, 

_ ^ec^,v-|z {Fu\z{u\z), Fv\z{v\z)\e) 
cu,v\z{Fu\ziu\z),Fv\z{v\z)\9) 

(11) 

Further, like for Ci^3|2, a C[/^\/|z term can depend on a parameter 9 through 
one of its arguments, say Fu\z{u\z, 9): 

d 

— In {cu,v\z{Fu\z{u\'^, ^yiz(^'lz))) = 

9cj7,y|z(J'c/|z("|z,6>),Fy|z(^^|z)) 

= — ( \ m F — ■ ^^t^|zHz>^ 12 

cc7,y[z (-?^c/|z(^^|z,6'),i^y|z(^^|z)) oO 

__ 9icuy\z (F^|z(^|z, g), Fv\ziv\z)) d_ 

" c^,v|z(i^^|zHz,0),Fv.|zHz)) -^/^l^^^l"'"^- 

Finally, in dimension c? > 4, both arguments of a Cuy\z copula term can 
depend on a parameter 9. In this case, 

— In {cu,v\z{Fu\z{u\z, 0), Fv\z{v\z, 9))) 

diCu,v\z{Fu\z{u\z,9),Fv\z{v\z,9)) d r i q\ 
- Cu,y,z{Fu\z{u\z,e),Fy,z{v\^,e)) • a^^^l-^^l^' ^) (13) 

d2Cuy\z {Fu\z{u\z,e),Fv\z{v\z,0)) d_ 
cu,v\z{Fuiz{u\z,e),Fy^ziv\z,e)) ' 89 ^l^^'^'^' >■ 

We see that the derivatives of copula terms corresponding to tree Tj in the 
vine will involve derivatives of conditional distribution functions which are de- 
termined by tree Tj_i. Thus, it will be convenient to store their derivatives in 
matrices 51*'""*'^ and SV'"^"'^^^'^ related to the matrices V'^*^"^* and y^'^^'-^ct 
which have been determined during the calculation of the log-likelihood to- 
gether with the terms 

In (cj{e)Me)\D{e){F3{e)\D{e){Xjl^e)\XD{e)),Fk{e)\D{e){^^ -■ Qj{e),k{e)\D{e), 
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for each edge e in the R-vine V , which can also be stored in a matrix V"™'"^'*: 

/... \ 



V 



values 



4,m6,5|'Tl7,5,m8,5 

Qa 

^?4,m8,5 ^3,m8,6 f?2,l 



(14) 



In particular, we will determine the following matrices: 



51 



direct.O 



-"41 



Qg F (u^lUjm , Urns, 5) gg F {u^lUrm ,s J Ums^e) 



^i^(u3|Ums,6) M^("2|mi) 



51 



indirect, 9 



(15) 



ae 



F{UI\U2) 



(16) 



51 



valuesfi 



ge f4,m6,5 my 5,7718,5 

9 ^ An 

m-7,5|m8,5 (96* ^^3, my, elms, 6 
9 „ d d 



(17) 



Here, the terms in Sl'^^^'^'^^'^ and 5'i«"-'^«''ect,e ^g^^^ determined by differenti- 
ating d?]) similarly as we did for the copula terms in (ITT]) - f ll3p . For instance, 
we have 



d_ 

de 



Fu\v,z{u\v,7,, 



d_ 



{hu\v,z {Fu\z{u\z,9),Fv\ziv\z))) 



d 



dihu\v,z {Fu\ziu\z, 9), Fv\ziv\z)) ■ — F[/|z(m|z,i 

d 

cu\v,z {Fu\z(,u\z,9),Fv\z{v\z)) ■ —Fu\z{u\z,9) 
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and 

d 



Q^hu\v,z{Fu\z{u\z),Fv\z{v\z,e)) = 

d 

= d2hu\v,z {Fu^z{u\z),Fv\z{v\z,9)) ■ —Fv\z{v\z,9). 

The complete calculations required to obtain the derivative of the log-likelihood 
with respect to one copula parameter 9 are performed in Algorithm 13.21 

Algorithm 3.2 Log-likelihood derivative with respect to the parameter 9^''^ . 
The input of the algorithm is a d- dimensional R-vine matrix M with 
maximum matrix M and parameter matrix 6, and a matrix C determined 
using Algorithm 13.11 for a parameter 9^'^ positioned at row k and i in the 
R-vine parameter matrix 6. Further, we assume the matrices V'^^'^'^'^^, 
ymdirect yvaiues corresponding to ouc obscrvatiou from the R-vine 
copula distribution, which have been determined during the calculation of 
the log-likelihood, to be given. The output will be the value of the first 
derivative of the copula log-likelihood for the given observation with respect 
to the parameter 6*^'*. 
1: Set Zi = v'^''^^"* 



k,i 

Set slfY"' := 0, s4"f := 0, sl^;^'"^^ := 0, t = 1, . . . ,d; k = i, . . . ,d 
if • == '■ then 



k,d—fh^ j+1 



Set Z2 = vf""""^ 

5: else 

6: Set Z2 = vf'^'''^'^ 

7: end if 

8: Set slf^'f* = dB--h{zx,Z2\B^~\9^~') 

9: Set sV^^lf = dB--h{z2,zi\B~^'\9~^'^) 

10: Set sUi'^''' - 



k,i expivVV^^") 

^ k.i ' 



11 
12 
13 
14 
15 
16 



for i = 2, . . . , 1 do 

for /c = A; — l,...,i-|-ldo 
if Cfc j == 1 then 

Q„f ^ ^.direct ~ „i direct 

sex Z\ — V^^^ , Zx — Sl^ • 

if mfc^j == mfc^j then 

Q„-j- ^ „, direct ~ „-\ direct 

oei Z2 — Wjfc Z2 — t'^k,d-rhk i+1 
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else 

c„f ^ „.indirect 5 „-\ indirect 

Z2 — Vk,d-rhk,i+l^ ^2 — ^i-k,d-rhk,i+l 

end if 

if Ck+i,i == 1 then 

Cpf „ -I values _ „-[ values 1 8ic(zi ,22 ~ 

iJtL ■'J-^^j — ^■^k,i ' "uoTues ■ ^1 

Set 51^;:^;=* = sifr.i';! + ^1/1(^1,^2 1 -b^'*, ^^'0 ■ 

Set sllT-iyr* = sl£lyr* + d2h{z2, zi\B'''\ 0^^') ■ zi 
end if 

if CkJri,d-rn+i 1 then 

Co+ eimities „-\ values i 820(21, 22|H'°'',g'°'') ~ 

•^e^ ^^k,i - ^^k,i ^^p^^values) Z2 

Set slf^.l';^ = slf:_fl + d2h{zi,Z2\B^^\ 9^^') ■ z^ 

Set sl^'L^^yf * = glindvrect ^ ^^^'^ " ^2 

end if 
end if 
end for 
end for 

return Ea^ ,=i d ^'^'I'f'"" 



In particular, this algorithm allows to replace finite-differences based numer- 
ical maximization of R-vine likelihood functions with maximization based 
on the analytical gradient. In a numerical comparison study across different 
R-vine models in 5-8 dimensions this resulted in a decrease in computation 
time by a factor of 4-8. 



4. Computation of the observed information 

Based on the calculation of the score function performed in the previous 
section, we will present an algorithm to determine the Hessian matrix corre- 
sponding to the R-vine log-likelihood function in this section. 
Again, we employ a convenient matrix notation. Considering a derivative 
with respect to bivariate copula parameters 9 and 7 associated with the 
vine, it is clear that the expressions for the derivatives of the log-densities in 
this case will contain second derivatives of the occurring h-functions. Thus, 
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our algorithm will determine the following matrices: 



52' 



\direct,d,f 



d 



-F{uA\Ums. 



■F{u^\ 



a) 



F{u,\ 



g<2'^ndirect,6,'y 
/... 



dea-y 



F(Ur 

d 



,M4) 



TOTS Pmg.s 



5, ''^4; 



96*97 



FiUmssWi) 



d 



F{u 



96*97 



F{u. 



'"■18,6 



(18) 



96*97 



(19) 



^'4,m6,5|m.7,5,m8,5 
'9 ^ 
9997'^4,m7,5|m8,5 
9 ^ 
96*97 t'4,m8,5 



9597 t'3,m.7,6|m8,i 
96*97 t'3,m8,6 



9 ^ 
909^ ^'2,1 



(20) 



Since not all entries in ([9]), (fTOl) and (fT4|) depend on both ^ and 7, not 
all entries in (fT8|) - (120|) will be non-zero and required in the algorithm. 
Employing Algorithm 13.11 to obtain matrices and C"' corresponding to 
the parameters 9 and 7, respectively, we see that the second derivatives of 
all elements where the corresponding matrix entry of either or C"' is zero 
clearly vanish. 

To derive an algorithm similar to Algorithm 13.21 which recursively determines 
all terms of the second derivatives of the log-likelihood with respect to pa- 
rameters 6', 7, we need to distinguish 7 basic cases of dependence on the two 
parameters which can occur for a term Cuy\'z (-F^|z(ii|z), -Fi/|z(f |z)) . 
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cases 


dependence on 


9 


7 


case 1 




Fu\z 


Fv\z 


case 2 




Fu\z, Fv\z 


Fv\z 


case 3 




Fu\z, Fv\z 


Fu\z, Fv\z 


case 4 


through 


Fu\z 


cu,v\z 


case 5 




Fu\z, Fv\z 


cu,v\z 


case 6 




Fu\z 


Fu\z 


case 7 




cuy\z 


cu,v\z 



Here, case 7 is relevant only for derivatives where 6 = '~f, since we assume that 
all bivariate copulas occurring in the vine density have one parameter in M. 
Because of symmetry in the parameters, all other possible combinations are 
already included in these cases, we only have to exchange 6 and 7. A more 



detailed description of the occurring derivatives is given in Appendix A 
As before, the terms in S2^''^'^^^'^'"' and S2'^^'^^^^^^'^''^ can be determined by 
differentiating (jTl) similarly as we did for the copula terms in ( lA.ll) - (1A.6[) of 



Appendix A Thus, the second derivatives can again be calculated recursively 
(see Algorithm Appendix A.l ). 

Combining Algorithm Appendix A.l and numerical integration technique^ 
we can also calculate the Fisher information (see Equation ([3])) matrix of 
R-vine copula models and determine asymptotical standard errors for ML 
estimates. 

Example 4-1 (3-dim. Gaussian and Student t-copula vine models) 

Let us consider a 3- dimensional vine copula model with Gaussian pair-copulas 
as given in the structure matrix M3 and the family matrix , ]Ye use the 

parameter matrix Q'-'"-^'^^ to denote the dependence parameters of the bivariate 
Gaussian copulas. 



'3 



B 



Gauss 



Gauss 
, Gauss 



e 



Gauss 



Gauss 



0.34 
.0.79 



0.35 



Then, we can calculate the asymptotic standard errors for each parameter 
based on the expected information matrix X(0'^""^*) or rather V'^'^p for the 



^We use the adaptive integration routines supplied by Steven G. Johnson and Bala- 
subramanian Na r asimh an i n the cubature p ackage available on GRAN which are based on 
Genz and MaUkI (|l980f) and lBerntsen et al.l ([199111 . 
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MLE rn fip (fipp. EqiiaU rm. | ! ^ ) ) and for the sequential estimation case (Equation 
and iHobcek Haft 1(201 1\) ), respectively. The order of the entries in the 
according asymptotic standard error matrices ASE^^^^ and ASE'^^'^ are 
the same as in the parameter matrix For 1(6^°'^'^^) and V^^p^ the 

order of the parameters is {pu, P23, Pi3\2) ■ 



j-M LE ^QGauss^ -jg^ 



^ QQGauss QQGauss 



Gauss \ 



Gauss\ 



'i,j = l,...,3 

0.65 
0.05 
-0.10 



1.62 -0.77 0.15^ 
-0.77 12.40 0.80 
0.15 0.80 1.42; 



0.05 
0.09 
-0.05 




This implies that the asymptotic standard errors are given by the square roots 
of the diagonal elements as 



ASE 



MLE 



0.86 

,0.29 0.80 



ASE 



seq 



0.89 

,0.31 0.83 



For the Gaussian distribution the occurring integrals can be computed analyt- 
ically, too, see\Appendix i?[ The results indicate that regardless of the applied 



estimation method approximately 100 observations are required to estimate 
the parameters of the 3-dimensional Gaussian copula up to a = 0.01. 
In the second setting we change the bivariate copula families to Student t- 
copulas. The corresponding copula parameters are stored m e^''"^''^\ where 
the lower triangle gives the correlation parameter (parameter 1 ) and the up- 
per triangle gives the degrees of freedom (parameter 2). As before ASE^^^^ 
and ASE^'^'^ denote the corresponding standard errors based on Equation 
Ij^ for the full ML estimation and Equation ^ for sequential estimation, 
respectively. 

( 3 3^ 

QStudent ^ Q 34 3 I ^21) 

\0.79 



0.35 



ASE 



MLE 



1.04 
,0.39 



12 



0.97 



ASE 



seq 



1.04 
,0.48 



12 



1.15 



14^ 
12 



Note that our results show that for uniform [0, 1] marginal distributions, the 
sequential procedure is less efficient than , full M LE. This is interesting in com- 
bination with Theorem 2 in lHobcek Han 1201 n) which states that together with 
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non-parametric estimation of the marginal distributions, the sequential esti- 
mation procedure for the Gaussian distribution is asymptotically as efficient 
as the full MLE. 



5. Application: Rolling window analysis of exchange rate data 



In this section, we apply the nietho ds de tailed above to the excha nge rate 
data analyzed by lCzado et al.l ( 120121 ) and lStober and Czadd (1201 if ). 
The data consists of 8 daily exchange rates quoted with respect to the US 
dollar during the period from July 22, 2005 to July 17, 2009, resulting in 
1007 data points in total. For simplicity, we use the following abbreviations: 
1=AUD (Australian dollar), 2=JPY (Japanese yen), 3=BRL (Brazilian real), 
4=CAD (Canadian dollar), 5=EUR (Euro), 6=CHF (Swiss frank), 7=1NR 
(Indian rupee) and 8=GBP (British pound). 



As ma rginal models we choose the time series models described by lSchepsmeier 
(l2010l . Chapter 5), which are of ARMA(P,Q)-GARCH(p,q) type. To ob- 
tain marginally uniformly distributed copula data on [0, 1]®, the resulting 
standardized r esiduals are tr a nsform ed using the non-parametric rank trans- 
formation (see Genest et al. ( 1995[ )). We could also employ the probability 
integral transforrn ation based on the parametric error distributions (IFM, 
Joe and Xu (Il996l )) but since we are only interested in dependence prop- 
erties here, we choose the non-parametric alternative which is more robust 
with respect to misspecification of marginal error distributions. We perform 
a two-step analysis, i.e. from now on we assume the marginal distributions 
to be known beforehand and only consider dependence analysis based on the 
obtained copula dat a. Given the marginal s, the sequential model selection 
procedure detailed in lDifimann et al.l (l201l[ ) suggests the R-vine described in 
the model matrix M (Equation ([8])) and pictured in Figure [1] as an adequate 
model for the whole dataset. More details on the involved bivariate cop- 
ula families can be found in [Appendix C[ The 8-dimensional R-vine model 



contains 28 pair copula densities; 8 of the conditional bivariate margins as- 
sociated with the vine can be modeled by independence copulas and for 7 
of the remaining margins two-parametric Student-t copulas are chosen (see 
i3. Equation ([231), with corresponding copula parameter estimates 0*^^^^ 
Equation (122]) ) Thus, the copula model has 27 parameters in total and 
we obtain an observed information matrix Xn{0^^^^) of dimension 27 x 27, 
where the parameters are ordered by column from lower right to upper left, 
i.e. Q^^^^ = (0.3,0.48,0.63,0.54,...). From this matrix, we can again ob- 
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tain standard errors SE^^^ (Equation (125]) ) for our parameter estimates by 
computing the square roots of the diagonal elements. As in Equation fl2ip . 
the estimates and corresponding standard errors for the degrees of freedom 
parameters of Student-t copulas can be found in upper triangular part of 
matrix (|22D and (125|). 



\MLE 



( 



11.85 





-0.70 


-0.09 




-0.88 


-1.44 


0.07 


-1.10 


-0.73 


0.26 


-1.23 


-1.17 


i^O.72 


0.55 


0.88 



1.13 
0.63 



1.08 
0.54 



0.63 
0.48 



0.3 



8.96\ 

7.74 

3.76 

9.97 

8.41 

7.46 

) 



(22) 



/ 

Indep. 
Indep. 
Indep. 
Indep. 

Gauss 
Student-t 
\Student-t 



Indep. 
Frank 
Frank 
r. Joe 
r. Gumbel 
Student-t 



Gauss 
Frank 
Frank 
r. Gumbel 
Student-t 



Indep. 

Indep. 
Gumbel 
Student-t 



Indep. 
Frank 
Student-t 



Frank 
Student-t 



Gauss J 
(23) 



MLE\ 



/1333.58 


-88 


15 


1 


39 


8.54 


-0.07 


-6.54 


-0 


25 


-88.15 


1843 


69 


5 


46 


-185.51 


8.41 


-27.61 


-1 


55 


-hi. 39 


5 


46 


27 


09 


-0.19 


0.04 


-0.28 





08 


-H8.54 


-185 


51 


-0 


19 


2279.43 


0.74 


-359.39 


200 


48 


-0.07 


8 


41 





04 


0.74 


29.21 


2.80 





75 


-6.54 


-27 


61 


-0 


28 


-359.39 


2.80 


3824.82 


374 


15 


-0.25 


-1 


55 





08 


200.48 


0.75 


374.15 


1803 


30 



7 

(24) 



20 



5.00 



2.54\ 
1.95 
0.59 
3.04 



SE: 



n 



MLE 



0.03 
0.03 
\0.01 



0.19 
0.19 
0.03 
0.03 
0.02 



0.03 
0.20 
0.19 



0.03 0.02 0.19 0.19 

0.01 0.02 0.02 0.02 0.03 



2.40 
2.01 



(25) 



We see that the parameters of selected Gaussian copulas are significantly 
non-zero. For all pair-copulas which are not the independence copula there 
is significant dependence and also the estimated degree of freedom parame- 
ters show that the corresponding conditional copulas are significantly non- 
Gaussian. To study possible parameter inhomogeneities in time, we apply 
a rolling window analysis as follows: A window with a window-size of 100, 
200 and 400 data points, respectively, is run over the data with step-size 
five, i.e. the window is moved by five trading days in each step. For each 
window dataset under investigation we estimate the R-vine parameters using 
ML while keeping the R-vine structure V given by ([8]) and the copula families 
B given by (123|) fixed. Additionally, we compute the observed information 
for each window and use it to obtain standard errors for the parameter es- 
timates. In some windows the estimated degrees of freedom parameters are 
very high leading to numerical instabilities and indicating that the Student-t 
copulas associated with some of the bivariate margins might not be appro- 
priate for the whole dataset. In Figure [2] we illustrate the estimated copula 
parameters together with pointwise approximate confidence intervals given 
by [9 — 2a,9 + 2a] for some selected pair-copulas. 

Our parameter estimates illustrate that the dependence structure of exchange 
rates in fact varied over the observation period although not as much as one 
might guess from a naive rolling window analysis without taking into account 
the estimation uncertainty. We see that the amplitude of observed variation 
depends on the window size with smaller windows leading to more promi- 
nent peaks. However, the greater uncertainty in parameter estimates due 
to less observations per sample leads to bigger confidence bands which will 
usually cover the parameter estimates obtained from bigger windows. Given 
our setting with a 8-dimensional financial dataset, 100-200 observations are 
insufficient to estimate copula parameters with satisfactory accuracy to de- 
tect short-term fiuctuations in dependence parameters (see Figure [21 row 5, 
with a peak in dependence visible in the results for 100 observations although 
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we have no indication that the parameter varied over time considering stan- 
dard errors and the results for larger rolling windows). Comparing results 
obtained with different window sizes, however, there is evidence that some 
dependence parameters vary over time. In the dataset at hand, the ob- 
servable changes in parameter values mostly occurred during the years 2007 
and 2008 in which the breakout of the subprime crisis lead to severe inter- 
ruptions in financial markets and changed paradigms in international cash 
flows. The increase of volatility is reflected in weaker dependence between 
currency pairs, e.g. the dependence between GBP/USD and EUR/USD 
decreased to 75% of its initial value in the year 2005 (see Figure [21 row 1). 
This is similar to the observations for stock markets, where dependencies dur- 
ing tim es of economic downturn are weaker than during times of economic 
upturn (IStober and Czadol (l201ll )). 



6. Discussion 



While ML inference for R-vine copula based models is broadly discussed 
in the literature and applied to data from various scientific disciplines, no 
algorithms for the computation of the score function or observed information 
have been available. 

The methodological contribution of this paper is to close this gap and to 
allow researchers to compute standard errors, which are at the core of ML 
analysis, in a routine manner for R-vine copulas. Furthermore, combining our 
newly developed algorithms with numerical integration techniques allows to 
compute the Fisher information matrix and other quantities related to the 
asymptotic theory of ML estimation techniques and thus to compare the 
efficiency of different estimation methods quantitatively for a given copula 
model. 

In a rolling window analysis we illustrate that the headline question of this pa- 
per is difficult to answer given the typical number of observations for datasets 
under investigation. Comparing estimation results for different rolling win- 
dow sizes and taking into account pointwise confidence intervals, we obtain 
clear indications that the strength of dependence between some currency 
pairs (e.g. GBP/USD - EUR/USD) varied over the observation period 
while it remained constant for others (e.g. EUR/USD and BRL/USD given 
AUD /USD). While the results for a window size of 100 trading days suggest 
that further short-term fluctuations in dependence might be present in ex- 
change rate data, these are non-significant in our modeling framework. For 
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the detection and study of such fluctuations, an analysis based on intraday 
data may be helpful and should be considered in future research. 
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Figure 2: Rolling window analysis for the exchange rate data with window size 
100 (left), 200 (middle) and 400 (right) for some selected par-copulas with 
different copula families (t-copula (row 1-4), Gumbel (row 5) and Frank (row 6)). 
The X-axis indicates the endpoint of each window, with the corresponding 
parameter estimate on the y-axis. The dashed horizontal line in each plot is the 
MLE corresponding to the whole dataset. 
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Appendix A. Algorithm for the calculation of second derivatives 

In Section m we introduced the seven possible cases of dependence which can 
occur during the calculation of the second log-likelihood derivative. In the 
following, we illustrate these cases in detail. In case 1 we determine 

ln(c,7,y|z {Fu\z{u\z,d),Fv\z{v\z,j))) 



ded-i 



- ' j^Fuiz{u\z,9)] ■ [—Fy\2{v\z,-f) 



cuy\ziFuiziu\z,e),Fy^^{v\z,j)) \de ^'"^ ' ' 'J ^dj 

- (^■^\n{cuy^2{Fuiz{u\z,e),Fy^2{v\z,^))) 

d \ 
{cu,v\z {Fu\z{u\z,9),Fv\z{v\z,'y))) \ , 



for case 2 



(A.l) 



ln(c(7,v|z {Fu\ziu\z,0),Fv\ziv\z,O,-/))) 



did2Cuy\z{Fuiz{u\z,e),Fv^z{v\z,9,^)) (d^ i \ A ( d ^ fin, 

- ' t:^Fuiz{u\z,9)] ■ [TrFv\z{v\z,e,-f) 



cu,v\z{Fu\z{u\z,e),Fv\z{v\z,e,-f)) \de ' J \dj 

d 

-^\rv{cuy\z{Fu\z{u\z,e),Fv\z{v\z,9,-i))) 
^H{cu.v\z{Fu\z{u\z,e),Fy\j,{v\'z,e,^))) 

+ / ' I I a ^\ ■ -^Fvizivz,e,-f) ■ —Fv\z{vz,9,-/) 

cu,v\z{Fu\z{u\z,9),Fviz{v\z,e,-f)) \d9 J \dj 

d2Cuy\z{Fu\ziu\z,0),Fy\ziv\z,e,^)) f ( \ o ) 

cu,v\z{Fuiz{u\z,0),Fy^^{v\z,e,j)) ^l^^''!^' 



(A.2) 
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and case 3 yields 



Hcu,v\z{Fu\z{u\z, e, 7), Fy\2(v\z, 0, 7))) 

cuy\z{Fu\z{u\z,9,-f),Fv\ziv\z,9,"f)) \d9 J \d-f ' J 

- In (c(7,y|z {Fu\z{u\z,9,-i),Fv\z{v\z,9,'^)))^ 

■ (^^H{cu,v\z{Fu\z{u\z,9,-i),Fv\z{v\z,9,^)))^ 

^ d\dicuy\z (f[/|z(M|z, 9, 7), j^y|z(f |z, 9, 7)) ^ ^ 
cc/,y|z (i^;7|z(w|z, 6', 7), Fv\x{v\z, 9, 7)) 



_^ - . - , j,v-|z (fc/|z(M|z, 6', 7), Fviz{v\z, 9, 7)) 

Ci7,V|Z (-Fi7|z(M|z, 9, 7), Fy|z(v|z, 61, 7)) 
_^ 9i92C[/,y|z (-F[/|z W|z, 61, 7), Fv\z{v\z, 9, 7)) 

cr/.v|z (f'r/|z(M!z, 9. 7), _Fy|z(i;|z, 6*, 7)) 

_^ 92g2Cj7,y|z (f j7|z(M|z, 6>, 7), j^y|z(f |z, 9, 7)) 

c;7,y|z (F,7|z(u|z,6',7),i^v|z(i'|z,6',7)) 

^ica,y|z {Fuiz{u\z,9,-/),Fviz{v\z,9,-f)) / 9 f I ~) 

cc/,y|z(i^c/|z(«|z,0,7),^^y|z(^^|z,0,7)) ' U^^T '^^ 

g2C[/,v|z (J'[/|z(H|z,g,7),-FF|z(^|z,g,7)) / 9' Mz9'y) 

Cu,v\z{Fu\z{u\z,e,^),Fy\z{v\z,9n)) \dedj ^'^^ ' ' 

Similarly, we have for case 4 that 



^Fuiz{u\z,9,j)^ (^-^Fuiz{u\z,e,'y)j 

a „ . , „ A / a „ . , „ A 



^i^y|z(t;|z,^,7) 

d 

— Fuiz{u\z,9,-f) 
d 

— Fv\z{v\z,9,-i) 

|2 



—Fu\z{u\z,9n) 



d 



—Fv\z{v\z,9,-i) 



d 



—Fv\z{v\z,e,^) 



(A.3) 



^2 

In (c[;,y|z {Fuiz{u\z, 9), Fviz{v\z)\j)) 
d 

^ln(c,7,v-|z {Fu\z{u\z,9),Fv\z{v\z)\^)) 



djdicu,v\z (f[/|z(M|z,6>),fy|z(f|z)|7) / / i 
cu,v\z{Fu\z{u\z,e),Fy\z{v\z)\^) ■ V 



^7Ct/,y|z (-F[/|z(M|z,6>),fy|z(i^|z)|7) 
cu,v\z {Fu\z{u\z,9), Fviz{v\z)\j) 



(A.4) 
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and 



^2 

In {cuy\x {Fu\x{u\7., 9), Fv\z{v\2^, 0)\l)) 



dedj 



^ In {cu.viz {Fuiziu\z, e),Fy\ziv\z, 0)h)) 

-djcuy\z (-Fj7|z(u|z, 0), Fviziv\z, 0)\-/) 
cu.v\z{Fu\ziu\z,0),Fv\ziv\z,9)\-/) 
d^dicuy\z {Fu\z{u\z,e),Fv\z{v\z,e)\-f) ( d 



(A.5) 



cc/,y|z {Fu\z{u\z,e),Fy\z{v\z,e)\-i) \d0 
d.yd2Cuy\z{Fu\z{u\z,e),Fv\z{v\z,0)\'-/) f d 



cuy\z {Fuiz{u\z,e),Fv\z{v\z,9)\^) \de 
for the fifth case. Finally, 
92 



Fu\ziu\z,e) 
Fv\ziv\z,e) 



In {cu,v\z {Fu\ziu\z, 9, 7), Fy|z(v|z))) 



didiCu,v\z{Fuiziu\z,e,j),Fy\ziv\z)) (d^ , ^ n .\ ( d ^ / i.n 

^ ^Fu\z{u\z,e,-f)\ ■ [—Fu\z{u\7,,e,-f) 



cuy\z{Fu\z{u\z,e,-i),Fv\z{v\z)) \de ^ ' ' y \d-i 
d 

In {cuy\z (-F(7|z(w|z, 7), Fv\z{v\z))) 

d 

In (cc/,y|z {Fu\z{u\z, 0, 7), i^v|z(t;|z))) 

^iCc/,v|z(j"a|z(^|z,g,7),i^v|z(Hz)) / 92 
c^,y|z(i^t.|z("|z,0,7),i^y|zHz)) ■^79^^!^^'''^''''^^ 
92 

^^ln(c,7,y|z (i^(7|z("|z),Fv|z(w|z)|6',7)) 

_ ded-fCuy^z {Fu\z{u),Fviz{v\z)\e,'-f) 
cuy\z (-F';7|z('"|z),i^v|z(w|z)|6i,7) 

_ decuy\z {Fu\z{u),Fy\ziv\z)\0,j) ■ djCuy\z (J^[/|z(m|z), Fy|z(^)|g, 7) 
cuy\z {Fmziu\z),Fv\z{v\z)\0,jf 



(A.6) 



Algorithm Appendix A.l Second derivative with respect to the parame- 

ters g^'' and g^-'. 

The input of the algorithm is a dimensional R-vine matrix M with 
maximum matrix M and parameter matrix 6, and matrices C'^'*, C^'^ 
determined using Algorithm 13.11 for parameters 6^'^ and of the R-vine 
parameter matrix. Further, we assume the matrices V'^"''^'^'^, and 
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■yvalues ^^-^q matrices j5']['^*''^'^*'^'* ^-^indirectjkji g^-Q^J ^-yvalueSjk,! g^-Qd ^\direct,k,i 

g-^indvrect,kj g-^vaiuesAi gj^^^^ rj.^^ output will be the value of the 

second derivative of the copula log-hkehhood for the given observation with 
respect to parameters 6^~- and Oj^-^. Without loss of generality, we assume 

that i > i, and k > k ii i = i. 



10 

11 

12 



if Cr - == 1 then 

k,i 

Set m = rhi 

Set zi = vfr--\ zi = slfr^*'^'^ 
if m —— i then 



17: 

18: 
19 
20: 



^ k,d-m+V ^ k,d-m+l 

else 

Set Z2 = v^f"^'^,, ~Z2 = sU^f ^'^^^•f'^ 

^ k,d-m+V ^ k,d~m+l 

end if 

Set s2f"'f? = 0, 521'^*';'='=* = 0, s2^'l^"^^ = 

fe— 1,1 re— 1,« re,i 

if A; == k i == i then 

Set s2f '■f * = (90- ,^0- M(zi, Z2\B~^'\ O'^'') 

ri — 1,2 k .i k .i 

Set s2|"_'';'|=^* = ,/i(z2, ^^'^ 



iq Clof cOvalues _ k,i k,i ^ ' ' ' ' _ values, k,i\2 

14: end if 

15: ifi'^. = lthen 

k+l,i 

Ifi- Cpi- „nvalues _ „'[values,k,i "k.i ^' ^' ' , "k.i ^' ^' ' ~ 



Set s2f ■^'^ = (9i(99- zzl-B'^'*, e^^) ■ z~i 

Set s2f_*Y'* = d2d0...h{z2, zi\B~^'\ 9~^'^) ■ zi 
end if 

if 4'\ , ==1 then 

k+l,d—m+l 



8280 c(zi,Z2\B 

Set §2^^ = g2^<^™^* -| 



Z2 



22: if c^:J.i . == then 

, i -80. -c(zi,22|e*'*,6>*'' 



fe,i fe,j exp(yl'4^'">^ 



24: end if 

25: Set s2f '^fl = s2f 7^ + 8280. Mzi, Z2\B~^'\ 9''^) ■ £2 
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26 

27 
28 
29 
30 
31 

32 
33 
34 



Set s2f^^f'* = s2i"_^ff * + dide- -.h(z2, zi\B~^'\ 9^^^) ■ £2 
end if 
end if 

for i = i, . . . , 1 do 

for k — k — l,...,i + ldo 
Set m — rhk,i 

if m —— ■mk,i then 

Oc>-t- y _ ^.direct zk,i _ -,direct,k,i ~k,i _ -,direct,k,i - _ 

oei Z2 — t^fc,d-m+l'^2 — *-^fc,(i-m+l' ^2 — ^2 — 

„rfdirect 
■^'^k,d-m+l 

35: else 

Qfi. Q„+ „ _ ^.indirect ~k,i _ indirect, k,i ~k,i _ indirect, k,i - _ 

6b. set Z2 - Vk,d-m+V^2 " ^^k,d-m+l '^2 " SJ-fc.d-m+l '^2 — 

„r)indirect 

end if 

fe,^ /Cji ' k — l]'^ ' ^ — 

if 1 & eg,,, == 1 then 

C!o+ c-O-^o-lues _ ^omZwes i didic{zi,Z2\B''-\e''-*) ~k,i ~k,i , dic{z-i,Z2\B'^'\9'^''') - 

^^k,i — ^'^k,i ^^p(_^values^ -Zl -Zi -f- ^ -Zl 

41: Set s2f:i'*^ = s2iZf;^ + dih{zi,Z2\B'''\9'''') ■ Zi + 

didih{zi,Z2\B'''\e''^')-~z'l^ -zl^ 
42: Set sT^^l^"* = s2*r*7^* + d2h{z2,zi\B^^\9'''') ■ zi + 

a2a2M;^2,;2i|i3^'\^''0-#^-5p 
43: end if 

44: if 1 & cgi_,_„+i = 1 then 

/It^. Q^-t- ^omiues _ „ovalues i a292c(2:i,Z2|g''-',6>'°-') , 92c(2i^2|^fe^ - 

45. DeX SZj.^. - SZ^^. -h • ^2 • ^2 + e^p(^.al„es) " ^2 



37 
38 

39 

40 



46 



Set s2fr^ = s2f:_\^l + d2h{zi,Z2\B'''\9'''') ■ Z2 + 

d2d2h{zi,Z2\B^^\9^^^)-~zl^ -zf 
47: Set s2f_*7^* = s2™_^^7^* + dih{z2,zi\B^^\9^^') ■ Z2 + 



didih{z2,z^\B'''\9^ 
48: end if 



~k,i ~k,i 
^2 ' ^2 



49: if c^ii, 1 & ct::M-n.+i 1 then 

en. Qo+ cO^'a'ues _ „^values , aia2c(2i ,22 ig'°'',6''°-') ~k,i 

ou. oei *Zj,^^ — AZ;^^^ -t- g^p(-^„atees) • ^1 • ^2 

51: Set s2f_fl = s2^!:^;f + did2h{zi, Z2\B^^\ 9^'') ■ zf ■ zf 
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52 
53 
54 



Set s2^"i{yf * = s2f:^;=* + did2h{z2, zi\B^'\ 6^^') ■ ■ 
end if 

if 1 ^ 4:;,. ==1 then 

56: Set s2f:^;=* = s2f;!^;=* + did2h{zi, Z2\B'''\ 6^^') ■ z^' ■ z''/ 

Set s2^"_^i7^* = s2f:f * + did2h{z2, zi\B'''\ 9'''') ■ z^' 



57 
58 
59 
60 
61 



end if 
end for 
end for 

return Efc,=i,...,. ^2-'«- 



Appendix B. Calculation of the covariance matrix in the Gaussian 
case 

While analytical results on t he Fisher information f o r the multivariate normal 



distribution are well known (IMardia and Marshall! (Il984j )) we will now illus- 





trate how the matrices Kq and (Equation and (^) can be calculated. 
We consider a 3-dimensional Gaussian distribution 

iV3(0,E), S= |pi2 

\Pl3 P23 

with density /123 and corresponding copula C123. Exampli gratia, we show 
the computation for the entry (2, 1) in Kg in detail. The other entries in Kg 
and ^g are obtained similarly. The first step is to calculate the following 
integral: 

■ln(ci2(ui,U2|pi2)) ) ( ln(c23(w2,W3|P23)) ) Ci23{ui, U2, U3)duidu2du3, 



'[o.i]3 \dpi2 J \dp23 

(B.l) 

where C12 and C23 are the corresponding copulas to the bivariate marginal 
distributions /12 and /23, respectively. Since the integral is independent of the 
univariate marginal distributions, we can compute it using standard normal 
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margins (see ISmithI (120071 )): 



12 



In(/i2(a;i,a;2|pi2)) 



dp: 



'23 



ln(/23(a;2,a;3|p23)) fi2?.{xi,X2,x^)dxidx2dxz, 



(B.2) 

where /12 and /23 are the according bivariate normal distributions. The 
3-dimensional and bivariate normal densities in ( 1B.ll) and ( IB. 21) can be ex- 
pressed as 



/i23(a;i,a;2,x3 



V2 



• exp 



' exp 



7r3/2^2pi3Pi2P23 ~ Pl3^ - Pl2^ + 1 - P23^ 

1 -Xi^ - X2^ - + a;i^/923^ + X2^Pl'i^ + X3^/Oi2^ + 2xiX2/0l2 +2x1X3/913 

2 -2 P13P12P23 + Pl3^ + Pl2^ - 1 + P23^ 
-2X2X3P23 - 2xiX2Pl3P23 - 2X1X3P12P23 - 2X2X3P13P12 



-2 P13P12P23 + Pl3^ + Pl2^ - 1 + P23^ 



and 



(B.3) 



fl2{Xl,X2) 



exp 



1 -Xi^ + 2X1X2P12 - 3:2" 



2^ Vl - I 2 (-1 + (Pi2 + 1 

Further, the derivatives needed in Equation (IB.2P are 
d 



(B.4) 



5p 

and 
(9 



ln(/i2(Xi,X2|pi2)) 



Pl2^ - a;iX2Pl2^ + - P12 + Xi^pi2 - X1X2 



12 



;-l + Pl2)'(Pl2 + l)' 



(B.5) 



dp 



Mf23{x2,X3\p23)) 



P23^ - a;2a;3P23^ + 3:3^23 - P23 + 2:2^23 - a;2a;3 



23 



:-l + P23)'(p23 + l)' 



Using ( !R3l) . ( !R5l) and (lR6l) in ( lR2l) we get 



(B.6) 



(US) 



Pi2^ - a;ia;2Pi2^ + a;2^Pi2 - P12 + 2^i^Pi2 - a:^ia;2 



(-1+P12)'(P12 + 1)' 
P23^ - a:2X3P23^ + X3^P23 - P23 + a;2^P23 - ^22^3 
(-1 + P23)^ (P23 + 1)^ 



/l23 (a;i , X2 , X3 )c?Xi C?X2C?X3 . 

(B.7) 
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The integral (1B.7P can be solved using well kn own results o n product mo- 
ments of multivariate normal distributions (see llsserlij (119181 )). 



(ElD 



+ 



(P23 + 1)'(-1+P23)'(-1+Pl22)' 

3 2 
— P12P23 + P23 P12 + P13 — P23 Pl3 



(B.8) 



(P23 + 1)'(-1+P23)'(-1 



P12' 



Since (p23 + 1)^ (—1 + P23)^ = (1 ~ P23)^ '^^'^ simplify Equation (IB.SP to 



(HH) = jESl) = 



(Pl3 - Pl2P23)(l - P?2)(l - P23) + 2pi2P23(Pl3 " Pl2P23)^ 



Pl3 ^ P12P23 
(1-P?2)(l-Pi3) " 
fcl2 

(1-P?2)(l-Pi3)' 



2/312P23 



(1-P?2)'(l-Pi3)' 
(Pl3 - Pl2P23)^ 



(l-p?2)'(l-pi3)2 



with 



= (P13 - P12P23) (1 + ^P^^P-^jr^^^0?^) ■ 



For the computation of terms corresponding to parameter pi3|2, note that 

Pi3 ^ P12P23 



Pl3|2 



and 



V(i-p?2)(i-py 



Pl3 = P13I2V (1 - P?2)(l - P23) + P12P23, 



which means that (IB. 31) . (IB.SP and (1B.6P can easily be re-parametrized. 
The final matrices are 



( 



fcl2 





fcl2 

(l-p^2)(l-p^3) 





\ 



1+Pl3|2 
W3|2-1)V 



\ 


1+Pl3|2 

\(P?2-1)(P?3|2-1) (P^3-1)(P^3|2-1) W3|2-1)V 



1+p!2 
(1-P^2)' 



Pl3|2Pl2 



23 

Tl2 



(i-pi3) 

Pl3|2P23 



35 



Appendix C. Bivariate copula densities 

This appendix introduces the bivariate copula families wh ich are included in 

our ex amples and applications. For more details we refer to lSchepsmeier and Stober 
fl2012h . 

The Gaussian copula with correlation parameter p G (—1,1) is defined by its 
density 

CGauss[Ui,U2\p) = ^^^-__ = exp | 2{1 - p^) J' ^ 

where Xi = 2=1,2 with $ being the standard normal cdf and 

(the quantile function) its functional inverse. Denoting the density of the 
univariate Student-t distribution with degrees of freedom h> as 



X 



2 



the bivariate Student-t copula's density is given by 

1 1 xl+xl- 2pxixi 

CStudent-t{ui,U2\p,iy) = == — — IH ^ 

277^/1 ~ at{xi,v)dt(x2,iy) \ uil ~ p^) 

Again, Xi is given as Xi := t~^{ui), i = 1,2, with t~^{-) now being the quantile 
function of the univariate Student-t distribution. 

In addition to these two well-known members of the elliptical class, we use 
the Gumbel copula, its rotated version to cover negative dependence, the 
Frank copula and the rotated Joe copula, which are so-called Archimedean 
copulas. 

The density of the Gumbel copula with parameter > 1 is given as 

CGumbeliui,U2\e) = C{ui,U2\e){uiU2)-'{{-Hui)f + i-Hu2)fr''^^ 

■ (In(ni) Hu2) f-' • {1 + (0 - 1)((- Hui) f + (- ln(^.2))')-^} 

with 

CGumbei{ui,U2\9) = exp[-{{-\n{ui)f + {~\n{u2)f}'^], 

being the corresponding cumulative distribution function. While the Gumbel 
copula can only cover positive dependence, it can be rotated to fit negatively 
dependent data: 

Crotated{Ui,U2) := c{ui, 1 - U2) . 
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Further, the Prank copula with parameter 6 e [—00, oo]\{0} has density 

CFrank{u-wu0) = ^(1 " e-'')e-^("i+"^) [(1 - e"^) - (1 - e-'^^){l - e-^-^)]-2, 
the density of the Joe copula with parameter ^ > 1 is 

Cjoe{ui, U2\e) = ((1 - Ui)' + (1 - U2f - (1 - - U2)') • (1 - Wl)'"^ 

• (1 - U2f-' . - 1 + (1 - Uif + (1 - U2)' - (1 - - U2f]. 
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